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Experiments on hybrid superconducting normal-metal 
structures have revealed that even in the absence of tunnel 
junctions the onset of superconductivity can lead to a decrease 
in the electrical conductance by an amount many orders of 
magnitude greater than e 2 /h. In this Letter we provide a the- 
ory of this phenomenon which shows that it originates from 
an instability in 4 - probe conductance measurements which 
is absent from 2-probe measurements. We compare the zero- 
bias, zero-temperature 4-probe conductances Gn and Gs of 
a normal diffusive metal in contact with a superconductor in 
both the normal (N) and superconducting (S) states respec- 
tively. In the absence of tunnel barriers, the ensemble average 
of the difference SG — Gs — Gn vanishes, in agreement with 
quasi-classical theory. However we also predict that there ex- 
ists macroscopic sample specific fluctuations in SG, which lie 
beyond quasi-classical theory and allow large negative values 
of 5G to occur. 



During the past few years, studies of the subgap con- 
ductance of hybrid normal-superconductor (N-S) struc- 
tures have led to the identification of a small number 
of paradigms of phase-coherent transport. The earliest 
of these is the zero bias anomaly observed in normal- 
insulator-superconductor (N-I-S) structures [ft) and the 
related finite bias anomaly found in N-S structures with 
clean interfaces. At high temperatures T > T* and 
bias- voltages V > V* , where for a N- metal of length L 
and diffusion coefficient D, k B T* = eV* = y^D/L 2 ), 
the conductances of both types of structure vary as 
\j\/T and 1/VV. For a clean interface there also ex- 
ists a conductance maximum at V* , T* and therefore at 
low-temperature and voltage a re-entrance to the low- 
conductance state occurs . This re-entrance phe- 
nomenon is also observed in a second paradigm of phase- 
coherent transport, which arises when a normal metal 
is in contact with two superconductors, whose order pa- 
rameter phase difference <fi can be varied by some exter- 
nal means. In this case the conductance is an oscillatory 
function of (p, with an amplitude which typically exhibits 
a maximum at T* and V* |2P-|l0[l . Unlike the Joseph- 
son current which decays exponentially with T/T*, such 
conductance oscillations decay only as a power-law. A 
third signature of phase-coherent transport is the long- 
range nature of these effects, which typically decay as 
a power-law in L/L*, where L* — y^D/eV). This be- 



haviour is in sharp contrast with the exponential decay of 
the Josephson effect and has been observed in a number 
of experiments [|ll],[l2| • 

All of the above phenomena can be explained by quasi- 
classical theory and to large extent are well-understood. 
There are however two further phenomena which are as 
yet unexplained by current theories of phase-coherent 
transport, namely the macroscopic suppression of elec- 
trical conductance by superconductivity |jb2| , [l3f and re- 
entrant paramagnetism in superconducting wires coated 
with normal metal O). Both of these effects have re- 
sisted detailed explanation, but must be understood if a 
complete theory of phase coherent transport is ever to be 
claimed. In this Letter we present a theory of the first 
of these phenomenon, which has been observed in several 
experiments |^,[b| , involving both non- magnetic (ie Sil- 
ver) and magnetic (ie Nickel) N-components, in contact 
with a superconductor via a clean interface. In contrast 
with quasi-classical theories |lf|[l£j which predict that 
when the superconductor is in its normal state, the zero 
temperature, zero-bias conductance Gat is identical to 
the conductance Gs in the superconducting state, these 
experiments show that the difference SG = Gn — Gs can 
be of order Gn and with a sample-specific sign. Since 
tunnel barriers and grain boundaries are absent from such 
samples, this effect is an embarassment to current theory. 
At finite bias, quasi-classical theory [|l6| has shown that a 
negative SG is possible, but the effect discussed in [16] is 
too small to account for the large changes found experi- 
mentally. In this Letter we suggest an explanation of this 
phenomenon which lies outside quasi-classical theory. 

To go beyond quasi-classical theory, we adopt a gen- 
eral scattering approach to dc transport, which was ini- 
tially developed to describe phase-coherent transport in 
dirty mesoscopic superconductors |L7| and is based on 
fundamental current- voltage relations derived in Jl8],[l9| . 
For simplicity in this Letter, we focus solely on the zero- 
voltage, zero-temperature conductance, for the structure 
shown in figure 1(a), which captures the essential physics 
of the experimental setup of reference fLl}] and comprises 
a disordered superconductor of condensate potential [i, 
sandwiched between two pieces of diffusive normal metal. 
The latter is in contact with four normal reservoirs at 
voltages Vj,j = 1, ... ,4. Leads 1 and 2 carry a current 
/, while leads 3 and 4 carry no current and therefore 
form the voltage probes of a 4-probe conductance mea- 
surement. As in the experiments, the superconductor is 
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chosen to be longer than the superconducting coherence 
length, so that sub-gap quasi-particles are not transmit- 
ted through the superconducting region. 

In the linear-response limit, at zero temperature, the 
conductance of a phase-coherent structure may be calcu- 
lated from the fundamental current voltage relationship 
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system and therefore in the presence of a magnetic field, 



v), 



(1) 



The above expression relates the current Ii from a 
normal reseviour i to the voltage differences (vj — v), 
where v = fi/e and the sum is over the 4 normal leads 
connected to the scattering region. The Q>ij s tire linear 
conbinations of the normal and Andreev scattering coef- 
ficients and in the absence of superconductivity satisfy 
Yl^j—i &ij — y^j—i a ij = in which case the left hand 
side of equation fil becomes independent of v. In units of 
2e 2 /h , an =Ni + Rf - R? and =Tfi-T%, 



where T? are Andreev and normal transmission co- 
efhccnts from probe j to probe i, Rf, Rf are Andreev 
and normal refelction coefficients from probe i and Ni is 
the number of open scattering channels in lead i. 

As an aid to comparing transport coefficients in the 
normal and superconducting states, it is convenient to 
eliminate v from equation ^ by the following procedure 
suggested in B3j2l| , which yields 
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Equation [2] is a fully quantum mechanical current - 
voltage relation which allows us to represent a multi- 
probe structure by an equivalent circuit, which for 4 
probes is shown in figure 1 (b). In contrast with the 
normal case however, the "conductances" Gij(i ^ j) are 
not necessarily positive. Furthermore, as noted in [ p0|j2l] 
Gij(H) — Gji(H*), where H is the Hamiltonian of the 



Gij 7^ Gji 



Setting I\ = —1% = J, J3 = J4 = and solving equation 
^ for the 4 probe conductance yields [Q 
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(G42G31 — G41G32) 



(7) 



where d > is the determinant of the 3x3 matrix ob- 
tained by removing the third row and column from the 
matrix Gy In what follows, we consider only the case 
where H is real (i.e. a real order parameter and no mag- 
netic field), so that G^ = Gji. 

As noted in [^TJ in the presence of disorder, the various 
transmission and reflection coeffcients can be computed 
by solving the Bogoliubov - de Gennes equation on a 
tight-binding lattice of sites, each labelled by an index i 
and possessing a particle (hole) degree of freedom ip(i) 
(cj>(i)). In the presence of local s-wave pairing described 
by a superconducting order parameter Aj , this takes the 
form 



Eip l = tiipi — Ytsl i^i+s + 4>i-s) + (Ai</>i) 
Efc = -ei4>i + J2 S 7 {4>%+s + <t>i-s) + A*Vi ; 



(8) 



The Hamiltonian of eq.(||) has been used to describe 
all of the paradigms [Til]. In the normal diffusive re- 
gion of figure 1(a), the on-site energy ej is chosen to be a 
random number, uniformly distributed over the interval 
eo — W/2 to e + W/2, whereas in the clean N-regions 

= eo- In the S-region, the order parameter is set to a 
constant, A$ = Aq, while in all other regions, Aj = 0. 
The nearest neighbour hopping element 7 merely fixes 
the energy scale (ie the band- width), whereas eo deter- 
mines the band- filling. In what follows we choose 7=1. 
By numerically solving for the scattering matrix of equa- 
tion 8, exact results for the dc conductance can be ob- 
tained and therefore if large negative values of 6G were 
not obtainable from such an approach, one would have to 
conclude that the experiments of [jl3| are not describable 
by mean- field BCS theory. In what follows, we employ 
an exact recursive Greens function technique to compute 
the scattering matrix of the structure of figure la. To 
obtain results which are typical of a normal metal, large- 
scale simulations are carried out, in which the sample 
dimensions are two orders of magnitude larger than the 
Fermi wavelength. 

Figure 2(a) shows results for ten different disorder re- 
alizations both in the normal state (Aq = 0) and in the 



As a simple example, for a symmetric structure where G13 
G24 G12 = G34 and G14 = G23 equation [?] reduces to 



2(G43 + G42XG43 + G41) 
G42 — G41 
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superconducting state (Ao = 0.1). As is seen from figure 
2(a), for certain disorder realizations the normal state 
conductance Gn is much greater than the conductance 
Gs in the superconducting state, yielding a large neg- 
ative value for SG = Gs — Gn- In contrast the en- 
semble average of SG, obtained by ensemble averaging 
over 500 disorder realisations, is negligible compared to 
Gn, in agreement with quasiclassical theory. As a pre- 
lude to the above simulations, to ensure that the disor- 
dered metal is diffusive, the conductance Gn was cal- 
culated as a function of length L with Ao set to zero. 
In the diffusive region, the mean free path I is indepen- 
dent of length and in units of the lattice spacing is given 
by I = (h/2e 2 )GNL/M, where M is the sample width. 
For the system used in figure 2(a), L = M = 200 sites, 
e = 0.2 W = 2.0 and I » 4.5. 

The fact that large negative values of SG are describ- 
able by equation ^ is a central result of this Letter. By 
examining individual scattering coefficients, we now ar- 
gue that this effect arises from the occurrence of small 
denominators in equation |?]. Indeed, in the normal state, 
if by chance G42G31 ~ G41G32 then Gn will be anoma- 
lously large. In this case, if the onset of superconductivity 
causes the denominator to increase, then a large negative 
SG can occur. 

Figure 2(c) shows the products G42G31 and G41G32 
for each of the normal state samples in figure 2(a). As 
an example, figure 2(a) shows that for sample no. 5, the 
normal state conductance is relatively large, while figure 
2(c) shows that G42G31 and G41G32 are approximately 
equal. Figure 2(b) shows the corresponding conductances 
when the samples are superconducting. As one might in- 
tuitively expect, the onset of superconductivity causes an 
increase in G42G31. Nevertheless the associated increase 
in the denominator of equation |?] leads to a large-scale 
decrease in the 4-probe conductance G. 

The aim of this Letter is to demonstrate that the phe- 
nomenon of large-scale superconductivity-induced con- 
ductance suppression in metallic mesoscopic structures 
can be understood through an exact solution of the Bo- 
goliubov - de Gennes equation. The results of figure 2 
demonstrate that for specific disorder realizations, as a 
consequence of an instability due to a small denomina- 
tor in equation [?], the onset of superconductivity can de- 
crease the total conductance, even though the individual 
conductances (Gy) may increase and even though the 
fluctuations in the Gy's are merely of order e 2 /h. Since 
the ensemble averaged total conductances for both the 
normal and superconducting state are equal, this phe- 
nomenon is not contained in a quasi-classical description. 

In addition to the calculations leading to figures 2, we 
have carried out extensive numerical simulations of two- 
probe structures and in no-case have we found large neg- 
ative conductance changes in the diffusive limit. Figure 
3 shows results for 10 disorder realizations obtained for a 
N-S structure of width M = 150 sites comprising a semi- 



infinte crystalline N-lead, in contact with a N-diffusivc 
region of length L = 30, which in turn makes contact via 
a clean interface to a semi-infinte crystalline S-lead. In 
all cases SG is negative and of order e 2 /h. A detailed 
analysis of the Lyapunov exponents shows that this is 
a consequence of the weak localisation effect predicted 
in p2j, which is too small to account for the large-scale 
effect described in this Letter. 

We believe that the results of figure 2 capture the 
essential physics of the experiments of Jl3| , although 
measured values for the conductances Gy are not avail- 
able, thereby making a quantitative comparison impos- 
sible. The crucial question is whether or not the condi- 
tion G42G31 « G41G32 is satisfied in the normal-state. 
Although the experimental samples are somewhat elon- 
gated, this condition may arise as a consequence of addi- 
tional disorder at the ends of the sample, where contact 
with the leads is made. For the future, it would be of 
interest to carry out an experiment on shorter samples, 
with a lcngth:width ratio less than or of order unity, since 
the denominator of equation [7] will be small. The sign of 
SG will be sample dependent and therefore by switch- 
ing on a magnetic field, fluctuations in Gy of order e 2 /h 
should generate both positive and negative macroscopic 
fluctuations in SG. Indeed there is no restriction on the 
sign of Gn and Gs and therefore such an experiment 
would also illustrate the prediction of superconductivity- 
induced negative 4-probe conductances discribed in ref- 
erence ptj. 
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FIG. 1. Fig 1(a) The N-S-N structur. 
equivalent quantum resistor network rcprcs( 
(a), where the Gjj arc given by equation 3. 
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2. 2(a) Results for the electrical conductance G N ,G g of the structure 
figure la, for 10 samples of differing disorder realizations. 2(b) Results 
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FIG. 3. Numerical results for 10 disorder realisations, obtained for the N-S 
geometry shown in the insert. 



4 



